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ABSTRACT 


Wet air oxidation (WAO) is one of the available technologies for the treatment 
of aqueous wastewaters containing organic pollutants. It enhances the contact 
between molecular oxygen and the organic matter to be oxidized. The reactor 
constructed on this principle contains deep well inside earth, which brings the 
reactants to the bottom of the reactor and takes out the by-products. The WAO 
requires certain temperature to start which makes the pre-heating of the reactor is 
necessary. 

This work is a validation of the formulation of the problem and its 
implementation as proposed in the earlier work by Mittal S. (1990). The earlier 
work presents a finite element model for the numerical simulation of start-up of 
the process inside a deep well reactor. It accounts for the heat transfer between 
the tubes and the heat losses to the surrounding earth and solves the heat 
conduction equation for earth and heat convection-diffusion for the fluid in the 
well tubes along with the mass balance of the fluid to find the time dependent 
temperature profile in the tubes and in the earth. Galerkin and Petrov-Galerkin 
variational formulation has been used for heat conduction equation of the earth 
and heat convection-diffusion equation of the well tubes, respectively. 
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Chapter 1 
Introduction 

Industrial revolution is a boon or a bane. This is a highly debatable topic. While it has 
brought upon a major change in human lifestyle in terms of comfort, it has also given birth to a 
problem that has put a question mark on the existence of mankind itself - pollution. Among 
the different kinds of pollutions, water pollution is a major one. It has been a very common 
practice of disposing the industrial wastewater into natural water bodies like rivers, lakes, 
ponds etc. This polluted water affects the population living around it both directly and 
indirectly. For this reason, treating wastewater before disposing it into natural water body is 
necessary. Proper and effective disposal of toxic industrial wastewater has always been a 
challenge for Environmental Engineers. 

Unfortunately, the options for effectively treating the toxic wastewaters and sludge are 
limited. Alternatives such as dilution, biological (secondary) treatment, physical-chemical 
(tertiary) polishing, deep-well injection, solar evaporation, incomplete or low-temperature 
combustion, and low temperature oxidation are not acceptable because the liquid and solid 
residuals may not meet future disposal criteria. Only two options - deep well-wet oxidation 
and incineration — are considered as effective as they are capable of destructing toxic organic 
wastes and sludge almost completely. 

Incineration with complete removal of stack discharges and ash stabilization has the 
potential of burning all organic material completely. But it has certain problems: (a) high 
construction and operating costs; (b) requires exceedingly high temperatures and complete 
mixing with assured detention; (c) generally limited to high concentrations of organic wastes 
in the fire box; (d) elaborate treatment trains involving considerable redundancy for treating 
stack gases, handling of ash and stack gas treatment streams; and (e) public resentment toward 
incineration of any kind. 
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The deep well-wet oxidation process has the potential to achieve many of the desired 
objectives. It utilizes the unusual property of water of being an excellent solver at high 
temperature and pressure. At high temperature and pressure, the solubility of inorganic salts 
is greatly reduced and the organic compounds are converted to harmless carbon dioxide and 
water. The kinetics of the process is controlled by reaction rates and not mass transfer. 
Relatively low waste concentrations can be treated due to efficient heat recovery. The 
environment is protected because the deep-well concept provides a totally enclosed treatment 
unit. Nothing is injected into the surrounding earth mass, left on the ground, or released into 
the atmosphere that is objectionable. 

1.1 Physical description of the system 

Deep well-wet oxidation is based on the use of a reactor vessel enclosed within drilled 
and cased wells. The design depth of the well may typically range from 1,524 m to 3,658m 
(5,000ft to 12,000ft). The reactor consists of a series of concentric tubes hung inside the earth 
from a surface platform. Fig 1.2 shows the configuration of a Deep Well Oxidation Reactor. 

The present work is a validation of the earlier work [1] done by Mittal S. So, the 
various configurations of the reactor and the parameters of study are kept same. In the original 
work, the reactor was assumed to have 4 concentric tubes. The detailed description of the 
reactor configuration is as follows. The concentric tubes, if we move from centre to radially 
outward direction, are known as steam tube, oxygen tube, down tube and the up tube. The 
steam tube brings hot steam or water to warm up the reactor. Down tube, as the name 
suggests takes the sludge down for the oxidation. Oxygen tube carries oxygen and the up tube 
takes all the residue of the reaction out. All the tubes end up in a common chamber called as 
reaction chamber. Reaction chamber is the place where the wastewater or sludge and oxygen 
get mixed and oxidation takes place. 



Fig 1.1: Schematic setup of a Deep Well Reactor 


The chamber is designed in such a way that the influents get sufficient time for getting 
oxidized under high temperature and pressure conditions. The outer most tube and the reactor 
are provided with multi-layered cement casing to prevent any possible seepage of toxic 
wastewater into the underground waterbed and thus provide extra safety from any possible 
contamination of underground usable natural water resource. 

There are many advantages of using deep well wet-oxidation. Since the water is 
introduced under a high pressure, the need for the high-pressure liquid pumps, heat 
exchangers and pressure-reducing methods that appear to be necessary for the aboveground 
systems, is eliminated. The operation of reactor is very economic as normally a single pass is 
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sufficient to provide 99.99 to 99.9999% removals of organic substances. Moreover it has 
sufficient flexibility that it can be operated at a wide range of operating conditions. If high 
destruction efficiencies are required, the reactor could be operated at low flow rates to provide 
sufficient reaction time for 99.99% destruction of organic compounds. If less destruction is 
necessary (80% COD removal), then the same reactor could process wastes at high volumetric 
flow rates. Lastly, since the reactor is totally enclosed, underground, compact and controlled; 
safety is inherent in its design. 

1.2 Applicability of the process 

The Deep Well Oxidation systems are applicable to a wide range of waste materials 
that include municipal primary, secondary, and mixed sludge; industrial sludge; high strength 
organic wastes; toxic liquid wastes; and solid hazardous wastes. These systems are designed to 
provide the desired temperature and density for the required reaction time. 

1.3 The Process 

The oxidation process cannot start at any temperature so the reactor must be preheated 
to the desired temperature. In the present model, high-pressure steam or water is injected 
through the steam tube for preheating. Hot oil can also be used for heating the reactor. 

Steam condenses due to heat exchange between steam tube and down tube and gets 
mixed with the down stream flow from the down tube in the reaction chamber. The oxygen 
supply is kept off during the start up process as no reaction takes place. 

Absence of any flow in oxygen tube results in a very small heat transfer coefficient 
between the steam and down tube. This prevents the steam from condensing entirely before 
reaching the bottom. The mixed stream goes out through the up tube. Similarly when the hot 
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water is used for start up, heat transfer takes place from down tube to up tube, which heats up 
the reactor tubes. 

Steam or water supply is switched off as soon as the initiation temperature of the 
reaction is achieved. In this model no extra pressure is applied as the hydrostatic pressure 
provides the sufficient pressure for the reaction. Once the desired temperature is achieved for 
start up, steam or hot water supply is switched off and oxygen supply is switched on for 
initiation of the reaction. Oxidation reaction, being exothermic in nature, fuels the further 
reaction between the sludge and the oxygen. Oxygen, at such high temperature and pressure 
conditions, mixes completely with water and organic materials, and with a near zero surface 
tension oxygen can penetrate minute pores. Thus, organic materials under such conditions can 
be oxidized to CO 2 and H 2 O very rapidly. 

1.4 The Finite Element Model 

For solving the problem numerically, several assumptions have been made: 

• Since the diameter of the tube is very small compared to the length of the tube, the 
temperature variation in radial and azimuthal direction has been neglected. So 
temperature remains uniform across a given cross - section and varies axially only. With 
this assumption, it is fair enough to assume the tubes to be one-dimensional. This 
automatically reduces the 3D problem to a 2D problem. 

• All the tubes are assumed to be frictionless. So the pressure in steam tube is assumed to 
be caused by hydrostatic head only. 

The model simultaneously solves the convection heat equations for all the tubes, conduction 
heat equation for the earth, mass balance equation of the fluid along with appropriate 
boundary conditions (discussed later) to find the temperature profile of the tubes and the 
earth. Since the fluid flow is one dimensional, Navier-Stokes equations are not solved for 
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determining the velocity field. The Petrov-Galerkin formulation has been used for spatial 
descretization of convection heat equation of the tubes to avoid any spurious node-to-node 
oscillations in the solution. The heat conduction equation for earth has been descretized 
spatially using regular Galerkin formulation. Predictor - Multicorrector time algorithm has 
been used to solve the set of linear algebraic equations obtained from finite element 
descretization. The details of the model are discussed in later chapters. 
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Chapter 2 

The Wet Air Oxidation Process 

In Wet Air Oxidation aqueous waste is oxidized in the liquid phase at high 
temperatures (400-573 K) and pressures (0.5-20 MPa) in the presence of an oxygen 
containing gas (usually air). The basic idea of the process is to enhance contact between 
molecular oxygen and the organic matter to be oxidized. High temperature conditions 
convert the organic matter to carbon dioxide and water. The liquid phase is maintained by 
high pressure, which also increases the concentration of dissolved oxygen (COD) and thus 
the oxidation rate. Typical conditions are 200-325 'C for temperatures, 50-175 bars for 
pressure and 1 h for the residence time; the preferred COD load ranges from 10 to 80 kg/ 
m3. 

2.1 Objective of the process 

It is an enclosed process, with a limited interaction with the environment. The 
process can treat any kind of organic aqueous waste that are too dilute to treat 
economically, even toxic wastes which can not be treated biologically, produced by 
various branches of industrial activity, or it can be coupled with a biological treatment 
facility to eliminate the sludge. WAO is one of the few processes that does not turn 
pollution from one form to another, but really make it disappear. This process has been 
used under mild conditions that alters molecular structures and improve biodegradability 
in spite of a limited COD reduction, but, usually, it is run under conditions that give 
almost complete oxidation of the organic matter and that can allow mechanical power 
generation because the reaction is exothermic. It has been established that the energy self- 
sufficiency is obtained provided that the COD load exceeds 12-15 kg/ m3. 



2.2 The Deep Well Wet Oxidation Method (DWWO) 

The deep well process employs wet oxidation method and is carried out in a 
vertical reactor buried in the ground. The reactor consists of concentric tubes of length up 
to 1500 meters or more suspended in a conventionally drilled and cased well. The waste 
enters the reactor at the surface and flows down the center tube,, which is referred to as the 
down tube or down comer. In the reaction zone, which is approximately the bottom 10-20 
per cent of the reactor, the waste stream and oxygen are brought together at high pressure 
and elevated temperature. Organic oxidation is initiated spontaneously at these 
conditions. Once initiated the exothermic reaction, combined with countercurrent heat 
exchange, fuels the process. The products leave through the outer annulus or up tube and 
exit the reactor at the surface. 

2.2.1 Purpose of the Well 

The purpose of the well is four-fold: 

1) To achieve the desired pressure in the reaction zone by the hydrostatic head. 

2) To allow large surface area for counter-current heat exchange between the 
descending reactants and ascending products so that the reactor can operate auto 
thermally (except during start-up and shut-down). 

3) To use the natural insulating properties and temperature gradients in the rock to 
minimize heat losses from the reactor. 

4) To minimize material requirements by taking advantage of the balance between the 
earth pressure outside and the reactor pressure in the interior. 
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2.2.2 Advantages of Deep Well configuration 

1) Safety: The pressurized heat exchanger/reactor is located below ground in a steel- 
cased well. Catastrophic failure of the reactor/heat exchanger will not result in a 
safety risk to operating personnel or other people in the surrounding area. No 
Carbon Monoxide or NOx is produced during oxidation. There is no threat of air and 
dust emission whatsoever because of its inherent closeness. The behavior of the 
reactor is always predictable because the mechanism of energy absorption is well 
known. 

2) Simplicity of design and construction: The reactor serves as the pressurization, heat 
exchanger and reaction vessel which eliminates recycling of supercritical water or 
use of high pressure, high temperature heat exchangers. The unit is constructed of 
threaded tubular pipes. 

3) High thermodynamic efficiency: The counter flow heat exchanger design enables the 
process to operate autogenically with a waste having a heat value as low as 1.2 x 10 
J/kg (50 Btu/lb). 

4) Broad operating range: The reactor is not confined to a single operating design point. 
Because of the large size of the reactor and the thermodynamic efficiency of the 
process, the unit can operate at different flow rates, waste strengths, and operating 
temperatures. 

5) Reduced pump horsepower: Hydrostatic head can produce the needed pressures. 

6) Earth serves as an insulator: Heat loss is relatively slow. 

7) Stability of operating conditions: The reactor volume is large. 
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8) Low operating cost: The processing cost for a deep well reactor in USA is roughly 10 
cents per gallon. This is the lowest operating cost among all the available options. 

2.2.3 Operational Flexibility of DWWO 

Despite being quite bulky. Deep Well Wet Oxidation process is flexible in operation. The 
process can be customized very easily, without much change in design of the reactor, for a 
given set of conditions. The process has following operational flexibilities: 

1) The pH condition in the reaction chamber can be changed from pH 2 to pH 8. 

2) The pH acidification can be by sulfuric acid or carbon dioxide gas mix in any 
proportion. 

3) The temperature in the reaction chamber can be set anywhere between 200F to 500 F. 

4) The reaction time can be readily changed from 2 seconds to five minutes. 

5) The flow rate can be varied 10 to 1 from a nominal rate of 5 ft per second. 

6) The concentration of suspended solids can vary from 0.01% to 14%. 

7) The reactor vessel can be cleaned while in service. 

2.3 Deep Well Wet Oxidation Process and Above the Ground Wet 
oxidation Process 

For both "top-of-the-ground" and "in-ground" processes, the following general principles 
of wet oxidation apply: 

1) The combustion reaction releases heat according to the following equations: 

C + a = ca + 14,100 BTU/lb. (from carbon) 

2 H 2 + O 2 = 2 H 2 O + 61,000 BTU/lb. (from hydrogen) 
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2) The combustion occurs in the aqueous phase. 

3) The heat of vaporization is conserved since the pressure in the vessel is always 
greater than the vapor pressure of the liquid waste at the oxidation temperature. 

4) The rate of oxidation is a function of temperature and characteristics of the waste. 

5) The degree of oxidation is a function of temperature and detention time. 

Apart from the principles stated above, the following process requirements must be 
satisfied for both "in-ground" and "top-of the-ground" systems, and must be cost effective: 

1) Presence of organic combustibles. 

2) Provision of high-pressure liquid, 

3) Provision of high-pressure source of oxygen; either pure oxygen or air, 

4) Effective heat exchange (for dilute i.e. less than 200 BTU/lb. waste), 

5) High-pressure reactor vessels with sufficient liquid detention time for the oxidation 
reaction to take place (20-30 minutes for sub critical, 1-5 min for supercritical). 

2.3.1 Advantages of Deep Well Wet Oxidation Method over above the ground methods 

Listed below are important advantages of deep well wet oxidation systems over 
equivalent units that are constructed on top of the ground: 

1) Use of Hydrostatic Head: The use of hydrostatic head of the reactor allows achieving 
the necessary reaction pressure by only circulating against the frictional head loss of 
the downdraft and updraft. This saves from 60-90% of the pumping cost, while 
significantly reducing the capital and maintenance cost of high-pressure pumps. 
This inherent process advantage also permits the deep well reactor configuration to 
process abrasive or other slurry wastes that cannot be reliably and economically 
pumped to extremely high pressures by top-of-the-ground system. 
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2) Combined Reactor and Heat Exchangers: Because of the concentric tubular deep well 
configuration the reactor acts as both the reaction vessel and the heat exchanger, 
thus eliminating the need for a separate external high-pressure tube-in-tube heat 
exchanger. The concentric tubular design of a typical 5000' reactor vessel provides 
excess heat exchanger capacity with nearly equal pressure on each side of the 
downdraft tube thereby reducing the cost of the downdraft/heat exchanger over top- 
of-the-ground external heat exchangers. For example, a separate heat exchanger 
would have to have 10,000 feet of 3" O.D. tubing of hastelloy material to be 
thermally equivalent to the deep well reactor "built-in" heat exchanger. Top-of-the- 
ground tube-in-tube heat exchangers encompass 180-degree bends every 40 feet. 
These are expensive and have been shown to be extreme points of erosion due to 
welded construction with un-annealed welds for the 180-degree bends. 

3) Natural Insulating Properties of the Earth and Use of Available Thermal Gradient: Thermal 
gradient provides a free 50-100 degrees F temperature increase for a typical 5000' 
reactor thereby allowing a more dilute waste to be treated. Additionally, the thermal 

• insulating properties of the earth provide an excellent natural insulator for the deep 
well reactor. The long-term storage of heat also reduces the heat required for startup 
after non-routine shutdowns, and provides potential for energy recovery after the 
useful life of the waste treatment project. 

4) Greater Thermal Efficiency: The greater thermal efficiency of the deep well reactor 
permits the treatment of a much more dilute organic stream than a top-of-the- 
ground system does without the use of a large heat exchanger or other methods to 
recover generated energy. Deep well reactor can operate with a minimum influent 
organic strength of 120 BTU/lb. compared to 350-450 BTU/lb. for top-of-the-ground 
sub critical reactors and 350 BTU/lb. for top of-the-ground supercritical reactors. The 
inherent greater thermal efficiency of these reactors, exhibit a much more 
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pronounced economic advantage for super critical application of low BTU wastes 
since the top-of-the-ground supercritical systems must employ recycling of 
supercritical fluid or employ expensive high pressure heat exchangers to recover 
heat. For example, a top-of-the-ground supercritical reactor treating a 10% organic 
(high strength, 1,750 BTU/lb.) waste requires a 50% recycle of supercritical fluid. This 
doubles the reactor vessel size for the same reactor detention time. For a dilute waste 
stream (2% organic with 350 BTU/lb. heating value), a 7,000-psi tube and shell heat 
exchanger of 1,000 to 1,800 SF is required for an 1,100 degree F reactor operating 
with supercritical water at a density of 10 lbs/CF. 

5) Fouling Potential: 

(a) Elimination of topside heat exchangers, which cannot operate in 2-phase mode 
(gas and liquid), is possible with temperatures (300 degrees F). This low 
temperature oxidation with sufficient oxygen availability, eliminates char 
formation and associated fouling from char formation in the reactor vessel/heat 
exchanger. 

(b) Comparatively low differences in temperature between reactor vessel influent 
and effluent means fouling from salts (calcium sulfate and calcium carbonate) 
that are inversely soluble with temperature are not as severe as with above 
ground units utilizing heat exchangers. 

6) Erosion/ Corrosion: 

a) Use of straight pipe with no bends eliminates erosion problems with U-bends in 
tube and shell heat exchangers. 

. b) Use of premium-threaded joints in the reactor vessel eliminated weld affected 
areas in reactor pipe that are susceptible to corrosion. All reactor pipes used are 
in fully annealed conditions to minimize stress corrosion cracking problems. The 
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bottom 100' plus or minus is made of extra thickness material to give extra 
erosion resistance. 

c) Lower reactor vessel heat exchanger temperature differentials reduce inorganic 
fouling over top-of-the-ground heat exchangers. 

7) Broad Operating Range of Reactor Vessels: 

a) Reactor vessels can be designed to operate over a 100% change in waste strength 
or flow rate to accommodate changes in reactor feed strength and volume. 

b) Low volume, high strength wastes can be treated successfully in reactor vessels. 

c) Due to continuous nature of reactor vessel/heat exchanger and long oxidation 
zone in the process, fouling effects, or change in oxidation rate of waste have less 
effect than with an above ground reactor with separate reaction vessel. For 
example, in a 5,000 feet reactor vessel, the upper 1/2 of the reactor, when it is 
clean, operates mainly as a heat exchanger, while the lower 1/2 operates as the 
reactor vessel. As the heat exchanger portion of the reactor vessel becomes 
fouled, the reaction zone moves downward into the reactor and, if the flow rate 
remains steady, the bottom hole temperature rises and the waste destruction 
efficiency remains complete. 

8) Safety: 

a) High-pressure reaction vessel/heat exchangers are located below ground in steel 
cased well. Catastrophic failure of reactor/heat exchanger will not result in safety 
risk to operating personnel or other people in surrounding area. 

b) If leakage from reactor or catastrophic failure from reactor occurs, all leakage is 
contained in cased well and is not release to the environment. 

c) If a dangerous situation arises, reactor and well can be flushed out and the 
material can be stored to prevent its immediate release to the environment. 
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9) Less Space Requirements: The deep well reactor vessel technology requires 
approximately 1.5 acres for an 80 DT/day reactor. This is approximately 1/5 the area 
required for commonly used alternative technologies. 

10) No Gaseous Emissions: The reactor vessels produce no odors or gaseous emissions 
that require EPA permits/ no NOx, no CO. 

2.3.2 Limitations of Deep Well Wet Oxidation 

1) It must have at least 25 tons of common organic materials to process during any 24- 
hour period, or 50 gallons per minute. It varies depending on the actual material. 

2) The reactor cannot be built within any strata that are subjected to displacement. 
These include salt domes and the crossing of a seismic fault line. 

3) The reactor is not portable in the sense that it is not truck mounted. It is relocate able 
and is designed for a twenty-year life. About fifteen percent of a typical 
capitalization cost would be lost to the abandoned cased housing chamber if no 
beneficial use for it can be defined. Suggested uses include mineral extraction and 
radioactive wastes internment. 

4) The reactor waste stream must not be so dilute compared to its volume as to fail to 
provide enough heat energy for both strata heat losses and the terminal temperature 
difference between the water entering the reactor and the water leaving the reactor. 
The lower limit appears to be 0.3 percent of organics in water with flow rates of 1500 
gallons per minute. Heat shortfall requires added energy. 

5) The reactor is limited by existing metallurgy to avoid peak operating temperatures 
in excess of 750°F if significant salts are present in the water stream. Even at cooler 
temperatures salt concentrations of 10,000 parts per million are suggested to be 
avoided or diluted. 
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6) When operating in an Oxidation mode the reactor feed should be limited to fewer 
than 9% of suspended materials in water. There are two considerations in this 
limitation. The first is the initial viscosity may be too great when dealing with sludge 
and result in too much pressure drop in the downdraft, which is counter to the 
preferred engineering. The second consideration is the release of too much heat 
energy in relation to the available water. Water is the principal heat-absorbing 
buffer. There is an operating mode where the excess heat may be converted to 
electric power or high-pressure heated process water but these are very special 
design cases. 


2.4. Types of Deep Well Wet Oxidation Method 

The deep well wet oxidation method is of two types: 

1) Sub critical deep well wet oxidation method 

2) Supercritical deep well wet oxidation method 

2.4.1 Comparison between Sub critical and Supercritical Deep-Well Reactor Designs 

The reactors based on the two processes have corresponding names. The design 
parameters and the optimum operating conditions for the two types of reactors differ 
significantly. Table 2.1 presents a comparison between basic design parameters of the two 
types of reactors. 



Sub critical Reactor 



'• - ; -it. % 

1200-1700 

2400-3658 


4000-5500 

8000-12000 
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' V; . A|S 

260-320 

400-510 

|f|g| 

500-600 

750-950 


//vppafc 

10000-14000 

22000-31000 

E V'.... .^r-Pv: P‘:£.U$ 

1500-2300 

3200-4500 

PH % W&xJ/v* ;'i.\ 

30-40 

0.1-2.0 


80-85 

100 

|v :;: r . ■■•'. ’"■ 

*- i. 'i. !»t S. AL,Juutiu LtoiJ&j! 

Oxygen gas 

Oxygen gas 

ppppsPy? 

Lower Vi 

Lower 1/20 

te : 

£ .' ^n.: .'liif .I.' 

1.2-27 

1. 5-4.6 


2-9 

5-15 

1 

Sludge, organic & toxic wastes 

All organic wastes 

— ~ — T-™ : 

»• .;• •’■'■•'r’f : 

Organic species dependent: 

80-99% 

99.99% 


Table 2.1: Comparison between Sub & Super Critical Reactors 
















Chapter 3 

Problem Statement 
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The aim of this study is to determine the temperature of the fluid in all the reactor 
tubes and the earth surrounding the reactor and below it during the start-up of the 
reactor. Further, we want to study the heat transfer from the tubes to the earth. This is 
important as the reactor may explode if the temperature of the reactor exceeds the safe 
limits. 

In this chapter all the governing equations and the boundary conditions needed to 
solve for the desired temperature profile in all the tubes and earth are discussed. 

3.1 Governing Equations 

The governing equations for the model are determined by the heat balance of the 
fluid in the reactor tubes and the earth and the mass balance of the fluid in the reactor 
tubes. 

Since the reactor tubes are cylindrical, the cylindrical coordinates would be the 
most natural coordinate system to work with. Given that the computational domain is 
symmetric in angular direction, it has been assumed that there is no angular variation in 
the temperature of the earth. Thus the standard Heat Conduction Equation in 2D with 
constant material properties gives the two dimensional heat balance equation for the 
temperature of earth: 

/ — \ dT}: 

^ X " p ) £ ~ir = 

Solving for the temperature of a flowing fluid usually requires the velocity of the 
fluid, which can be found by solving the Navier-Stokes equation for fluid flow. We have 


l d( dT£] d% 
r dr V dr J dz 2 


( 1 ) 
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assumed the tubes to be one-dimensional and the fluid velocity to be uniform in all the 
tubes. So Navier Stokes equation has not been solved. 

The 1-D schematic of the tubes is shown in fig 3.1 



Figure 3.1 Schematic of the reactor 


The heat balance in the steam tube results in 


l A PCp%f§- + 


pm 


dx 


m s php 


+ (ouy d {r§-T$)=o, o<x<h 


( 2 ) 


V J \ J 

Since we have assumed no friction inside the tubes, the pressure distribution would be 
hydrostatic. This is given by 


dpp s 

fp = PF8' 0<:X<H 

The mass balance in the steam tube gives 


(3) 


b dt dx 


0<x<H 


(4) 
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The heat balance equation for the down and up tubes are given by 


[ApC f Y+{ApC P %0- + [mCp) d ^--(aUr^-T^-faUr^-T^S, 0<xiH 


(ApCp)p + (ApCp)p- 


+ m 


mCp) U ^ + (oUY d {n -T$) = Q, H<x< 2 H (6) 


To ensure the mass continuity at the reactor, the following equation is used. 

• • • 

m.p(t) = m d p + m S p(H,t) (7) 

This equation is used to determine the mass flow rate in the up tube. The mass flow rate 
of the down tube is considered as a constant. 


3.2 Thermodynamic Relations 

The governing equations described above cannot be solved completely (with 
given boundary conditions) because the enthalpy and density terms in the equations are 
functions of temperature and/or pressure. So we need the proper thermodynamic 
relations. For single phase the enthalpy and density are function of both temperature and 
pressure. 



Using the definition of complete derivative, equation 8(a) can be differentiated with 
respect to time and the ID space coordinate to get 
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dh b F 

dx 


<%} 

V . dT Fj 


Vf 


El 

dx 

\ J 


+ 


dhp 

dp s F 


t: 


dx 

\ J 


9(b) 


Now using 9(a), 9(b) and (4), the heat balance equation for steam tube for a single-phase 
fluid can be modified to 


(ApCpl + A s F p S p 


8 hr 


dTv 


'V r 


dTp 

dt 


• + trip 


Shi 


dTr 


dT F 

dx 


' (oiiyd^Tp — = -I 


' V, 


dtip 
dp S F 

v j T s 


/a \ s ^Pf s d Pp 


0<x<H 


( 10 ) 

For a 2-phase system (saturated water), temperature is considered as a function of 
pressure 



And the quality of the steam is used as an independent variable. The specific enthalpy of 
the water steam two-phase system can be expressed as a linear combination of specific 
enthalpy of water and specific enthalpy of steam. The weight function depends upon the 
quality of the steam. This relationship is given by 

4 = 44 + ^- 4 ) ( 12 ) 

Moreover, the specific enthalpy of steam and water are functions of temperature 
themselves 

K=K{ t f) 13(a) 

h i =h i{jF^ 13(b) 

So equation (12) gives the specific enthalpy of the steam at a given temperature. Similarly 
the density of the steam can be determined by following relationships for a given 
temperature of steam 


vrrs-jArr 
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Pp - p v Sp + Pi l-£p 


Pv Pv , 


15(a) 


Pi * Pi[ t f 


15 (b) 


Since the enthalpy is a function of temperature and quality, the following relations hold 
good 

I'M) 


dt l dTp A dt 


dhp _ dhp 


& [dTr I . dx 


16(b) 


Now by substituting Eq 12 in 2 and using 4, 16(a) & 16(b), the heat balance equation for a 
two phase system in the steam tube can be written as 


(ApC P Y + A s f p s f \ —L. 


dTr % dhi 
— - + mr — — 

dt dTp 


+ (oU) sd lTp-Tp\ = 


dht s del % del 

< )r c 


where 


F dTp ■ Arl 


0<X<,H 


A A 

K~ h i 
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3.3 Initial and Boundary Conditions 


The partial differential equations are unsolvable without boundary 
conditions. The various initial and the boundary conditions used for the given model are 

described below 


Initial Conditions 


Tp(x / 0) = T E (x) 0<x<H 

(20) 

Tf(x,0) = Tp(x) 0<x<H 

(21) 

Tp(x,0) = T°(x) 0<x<H 

(22) 

Tp (x,0) = Tg (2H - x) H<x< 2H 

(23) 

Boundary Conditions 


T E (r 0 ,z,t) = Tp(2H -z,t) 0 < z < H, t > 0 

(24) 

T E (r 0 ,z,t)=Tp(2H-z,t) 0<Lz<H + h, t>0 

(25) 

T E (r,H + h,t) = T E (H + h ) 0<r<h, t>0 

(26) 

3T 

— —(0,z,t) = 0 H<z<H + h, t>0 

dr 

(27) 

^ {r o ,0,t)=0 r 0 <,r£h,t>0 

(28) 

Tp (0,t) = T E (0) t>0 

(29) 

Tp(0,t) = Tq( 0) t>0 

(30) 

T E (r / z,0) = T E (z) 0<z<H, r Q <r<h and H <z<H + h, 0 <r <h, t>0 

(31) 

• • 

m S p(0,t) = rn s 0 (t) t>0 

(32) 

• m 

md p{O r t) = m d 0 (t) t>0 

(33) 

Tp(H,t) = Tp(H,t) = Tp(H,t) t>0 

(34) 
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p s F (0,t) = pQ(t) t>o 

£ S F {0,t) = 4( t ) t ^ 0 


(35) 

(36) 


3.4 Non-Dimensionalization 

For Non-Dimensionalization of the governing equations, following scales have 
been used 

T max- T E (0) (AT) max 
* _ T * _ X 

r ~h' x ~H 


Z = 


H 


1 + 


z-H 


0<z<H 
H<z<H + h 


t = 


a = 


JpCp)r 


The governing equations in their Non-Dimensional form are given by 


80 _ 1 8 

_ * * ^ * 
dt r dr 


80 


r 

\ dr" 


1 d~0 E 

7^ 


+ (Vf (4 -ety (Vf (4 -4) = Z" ° 

dt 1 8x vyv J J \ r t ) 


< x <1 


(37) 


(38) 


dt ij ctc v/v ' '• 


0< x <1 


1< x <2 (39) 


(40) 
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where, 

( A pc P Y p 

(AfiCpi 




{ApC p Y p + {ApC p Y F 

WCpf F 






for single phase 


( 41 ) 

( 42 ) 

( 43 ) 

( 44 ) 

( 45 ) 

( 46 ) 

( 47 ) 


for two phase 


( 48 ) 



(49) 



h 2 {aU) ud 

cc{ApCpYp 



>i(ouy^ 

a(ApC p Y F 


h 1 2 S 

a (^T) max [ApCp Y F 


S = 


S S , A A 

m F £ F h v -h t 


D 




(50) 

(51) 


(52) 


[ApCpYp 

hjc 

a (ApC P Y F 


(53) 

(54) 



1 = 1 z*>l 

L=*L OZz*£l 0< x* £2 
h 


for single phase 

(55) 

for two phase 

(56) 


(57) 




(59) 


(60) 
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for single phase 


for two phase 


(61) 


(62) 



Chapter 4 

The Numerical Solution 


28 


The governing equations described in the previous chapter are non-linear partial 
differential equations. Generally, it is not possible to find a closed form solution for non- 
linear PDEs with given boundary conditions. So the governing equations are solved 
numerically with the boundary and initial conditions described in chapter 3. 

Algebraic equations of the form Ax=B are far easier to solve compared to PDEs. So 
the Numerical methods are used to reduce PDEs to algebraic equations. In this study we 
have used FEM for this purpose. The algebraic equations thus obtained are then solved to 
obtain the desired solution. 

4.1 Spatial Descretization 

The whole computational domain has been divided into four computational sub- 
domains; one for earth and rest for the reactor tubes. The computational domain for earth 
Qe can be descretized into sub-domains called elements given by D% e - 1 , 2, 3 .... nek. 
The number of elements, nek depends completely upon the physical set up of the problem 
and the desired accuracy of the solution. In the present study, the computational domain 
is rectangular in shape and the elements are also rectangular. The computational domains 

for the reactor tubes £2p,f2p, f2p are linear. So the elements of these domains are linear. 
Because of the difference in dimensions of the elements of earth computational domain 
and the reactor tubes computational domain, the implementation of various FEM 
techniques is very difficult. So for the reactor tubes computational domains, two- 
dimensional pseudo-elements have been used. Each pseudo-element consists of two one- 
dimensional elements facing each other. This has been depicted in figure 4.1 

a AL 
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Figure 4.1 Two-dimensional Pseudo element 


X 

4.2 Variational Formulation 


For setting up the variational formulation we need to define sets of weight 
functions and the solution functions. These finite dimensional spaces, can be defined as 
follows on the basis of the boundary conditions. 

V h (n E )= w h w h e H lh (l2 E ),w h = 0, onT (63) 

o 



(64) 


The two spaces can be called as Variation Space and Solution space respectively. 

The discrete variational formulation (Galerkin Formulation) associated with the 
heat conduction equation of earth can be given by the following equation 


r n UUr * * * 

jzu — &-r dr Ldz + J 


dw dd E 2 dw h ddr i * * * 

L 1 r dr Ldz 


oAdr dr 


L d z dz 


Flere we wish to find all d h £ S h by solving above equation such that w A £ Y* 
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The heat flux at the up tube - Earth interface is approximated over a pseudo-element 
(shown in figure 42) as 


Irry 


l dr* 


r =r n 


= 2jry 


* Ar 

Tq h 


r 




6_-6 
v J 


( 66 ) 



The computational domains for the 


Dp = \x 


0< x <1 


I2p — I x 


0< x <,1 


Dp = I x 


1< x <2 


reactor tubes or fluid can be defined as follows 

67(a) 

67(b) 

67(c) 


The Variation and Solution spaces for these computational domains can be defined as 
follows 
Down Tube 

v ih [4) = 


dh 


w 


dh 


w 


e H lh \n 


ZD 


dll 


0 at x = 0 


68(a) 
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s ih ( ai) = 

Steam Tube 



v sh (n F 

s sh [n s F ^j 

Up Tube 




| e F dh e H lh (n F y 0 F dh = 0 at x* =0 

w sh e H lh j, w sh = 0 at x = 0 
e F sh g H lh [n F y e F sh = e s 0 (t) at x* =o 


68(b) 


69(a) 

69(b) 


v uh (n u F | = 


uh ( 


uh 


TV 


uh , r lh 
\W gH 


(^ 


12 




GH“|fl 


lh{ ^u 
F 


70(a) 

70(b) 


To ensure proper heat transfer at the bottom of the reactor the Solution Space has to 
satisfy following constraint 

(71) 


ef =ef =9? at x* = i,\tt>o 


The Galerkin discrete variational form is known to generate spurious node-to- 
node oscillations for convective equations. So for the convective heat balance equations 
for the reactor tubes Petrov-Galerkin method has been used. The variational forms 
associated with these equations (38, 39 and 40) are given by 
Down Tube 
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Up Tube 


n uh c.uh 

w + 8 
nl 


k} 


def {Pef ddf 


dt 


rMh ~dh 




+ 2ny 


( *M 
* Ar 




(73) 


dn=o 


Steam Tube 

.( sh ,sh4l f \d0p {Pef ddf ( ,\s d( sh dh \ 

J w + <? }Vpf* + W - + {P \ d F ~ 9 F } dn 

n $ p \ dt l gx \ j \ ' 

= j (w sh + s sh ^p*) s dn 


cr r 


(74) 


Where &*• is the Petrov-Galerkin supplement function and it is defined as 


, .h ~ Ax 

5 = C 2 r— St 8 n \ 


( 71 ~ \ Zu.J 1 - 


Pe 

\fpj 


dw 


dx 


Where C 2r is the Courant number. For this study, Courant number is unity 


C 2r=I 


(75) 


(76) 


Using the value of Courant number and writing the differential in the local coordinate 
reduces the supplement function to 


8 = signl 


\fpj 


dt 


(77) 


Let c= sign 

So, finally, 

~h dw h 
8 =c- 


r Pe^ 


\fp J 


H 


(78) 
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4.3 Elemental Matrices 


The variational formulation of the governing equations of the reactor tubes gives 
rise to the elemental matrices. The elemental heat capacity matrices, convection matrices 
and matrices corresponding to the heat exchange terms are given by 



( ■' \ sd 
X, 


Ax 




^2 c ^ 


(2 C 1 1 

f 1 


f 1 c N 

U 2, 


1.3 2 ) ' 

^3 2 j 


l3 _ 2y 


(1 c 

1 3 2y 


(1 c 

j 

U 2 


2 c 

1 

U 2 


fi+i) 1 

U 2 J 1 

r 2 c ^ 
L 3 + 2 y 


f 1 c ) 
L 3 + 2 J 


^ 2 c 'j 

ri c 


(2 cS 


ll'lj ' 

U 2J 


U 2 J 



f 1+£ 

U 2 
r 2__c_ 

{3 2 

(2 c 

U 2 


(81) 


The elemental force vector due to heat source in up and down tubes are given by 


(/ f f =4 


^b,* + 

3 2 


.3 2 


2 c\ * f 1 c \ * 

+~ IX + V+-T 


3 2 


3 2 


3 2) 2 U 2 ‘ 

(1 c ^ „* (2 


U 2 


~+~ x + 


“I x 

J Z 


(82) 
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The elemental force vector for the steam tube is given by 


A 


sd Ax* 


---Yp*Y+(---Yp*T 

3 2 A h U 2 A k 

0 

1 + iYp-Y+ri £ v .y 

3 2 A h U 2 A k 


(83) 


4.4 Assembly 

The assembly procedure is described by the following equations 


p _ A p e , A 
L ~eeE £ L E + eeEp 


Y — A y e . A 
K eeEpE + ee£p 


r A f e , A 
t eeEpJE + eeE F 


A 


FJ ~ r \FJ 


+ c 


A s 


4 


\ud 

■■ \ud 

( 

sd 

•• \ sd 

+ 

4 

+ 

4 

+ 

4) 

) 

i ) 


V J 


i j 


JL. A V e 
+ eeE ^ F 


w 




sd 


(84) 

(85) 

( 86 ) 


4.5 Solution of the simplified equations 

The assembly of the elemental matrices leads to a set of time-dependent ordinary 
differential equations. 

C — + Kd = F 

dt (87) 

d(0) = 0 

Or, 

Cv + Kd = F 88(a) 

d (0) = 0 88(b) 

This system of equations has been solved using Predictor-Multicorrector method. The 
solution procedure has been described below 
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The solution for n+lth step is predicted by 
Pn+l =d n + t 1 " 1 ) t V n 


And corrected by 

d n+l = V n +l +;itv n+l 
Where 


(89) 

(90) 


v n+l {C n+ l^tK n+i y [f n+ i K n+iPn+l ) 


(91) 
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Chapter 5 

Results and Discussion 


The start-up of the reactor is studied for the case when hot water is introduced in 
the steam tube. The following numerical quantities were used to set up the physical 
problem 
ro = 0.52 '.ft 
H = 5000 ft 
h = 10ft 

Tg(z) = 80 + 0.628z °F 
(pC p ) E = 40 BTU/ft 3 
(a) p = 0.037601277 ft 2 
(A); = 0.097555518 ft 2 
(A); = 0.012271846 ft 2 
(A) P = 0.132822283 ft 2 
(a) u f = 0.20141904 ft 2 
= 0.021816615 ft 2 

{ApC p ) d F = 8.234981546 BTU/ft °F 

a ud = 1.65915362ft 2 / ft 

a sd = 0.916297857 ft 2 / ft 

Pj = 1201.5341 

PJ = 2.748893571 

f =71.8614893 

//= 0.160667625 



//= 0.416847903 
// =0.052436739 
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A mesh of 2813 nodes, 2536 two-dimensional elements and 228 one dimensional 
elements (152 pseudo-elements) has been used to descretize the computational domain 
which is a rectangle. The total number of 2688 two-dimensional elements has been used. 



c 



\ 



The reactor was started up by introducing hot water at 500 °F and 700 bars 
through the steam tube. Actually the temperature at the inlet of the steam tube has been 
ramped up from 50 °F to 500 °F in 1.2 hours. The mass flow rate at the inlet of steam tube 
is assumed to be 100,000 Kg/s and inlet of down tube as 50,000 Kg/s. Figure 5.1 shows the 
temperature profile of the tubes after 12 hours of introduction of water. This temperature 
profile is in agreement with the previous work [1). 
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The bottom temperature (bottom of the reactor tubes) is a parameter to determine 
the time to start up the reaction. Hot water continues through steam tube for that time 
only for which the bottom temperature has not reached the desired initiation temperature 
of the reaction. Once the initiation temperature is reached, hot water supply is shut off 
and oxygen supply is started. The bottom temperature variation for 12 hours is shown in 
figure 5.2. The bottom temperature first decreases rapidly, then increases sharply and 
finally becomes almost constant. The final bottom temperature is very close to the inlet 
temperature of the steam tube. 

Figure 5.2 

Bottom Temperature Vs Time 


The fall and rise of the bottom temperature has been captured in figure 5.3 


Figure 53 Bottom Temperature Vs Time 
Initial part of the curve 



Time (Hours) 
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Again the variation of the temperature is same as found by Mittal S. (1990). Another 
important parameter to study is the temperature of the fluid at the exit of up tube. It is a 
measure of actual heat transfer from fluid to earth. The exit temperature increases 
constantly with the time as heat is being added to the system through the inlet of the 
steam tube constantly. After a long time, at steady state, it becomes flat. Figure 5.4 shows 
how the exit temperature varies with time. 



Figure 5.4 Exit Temperature Vs Time 


The exit temperature variation is also found to be the same as the previous study [lj. 




Chapter 6 

Numerical aspects of the Wet-Oxidation Reaction 
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Once the reactor is heated to the desired value of bottom temperature, the hot 
water/steam supply is stopped and the oxygen supply is opened. The oxidation reaction 
of the sludge is exothermic in nature; the heat produced by the reaction carries on the 
process further. In this chapter, a possible method for the numerical study of the reaction 
has been suggested. 

6.1 Governing Equations 

The governing equations of the fluid described in chapter 3 still apply but now we 
will have conservation of species law in addition to the conservation of mass, momentum 
and energy. 

• • 

Let m d be the total mass flow rate in the down tube. The mass flow rate m i and 

mass fraction co i of the i* species are related as 

= ™d a i 

Let p be die local fluid density, p { and r- be the density and rate of production of i* 
species. So, 

Pi = P°>i 

If Ad is the area of cross section of the down tube then the continuity equation for the i* 
species is given by 
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The Total Continuity equation will be the sum of the continuity equations of all the 
species. Since = 1, the Total Continuity Equation then becomes 


A d — + = 0 

1 dt dz 


Since very dilute sludge is used in Deep Well Wet Oxidation Process, it is 
reasonable to assume that the local density of the fluid does not vary much with respect to 

time i.e. ^ » 0 . The Total Continuity Equation then becomes 


dm l 

dz 


■■0 


i.e. the axial variation of the mass flow rate of the i* species is negligible. In fact the axial 
variation of mass flow rates of all the species would be negligible. If there are N-l species 
then for the i* species the continuity equation becomes 


dco: m„ dco- 


dt A d dz 


l _ 


= 'i ,i=l2 N-l 


At the bottom of the reactor continuity is assured by 


m u= m d+ m o , 


The continuity equation in the up tube is given by 


da); m„ dco- 

P ~dt + ~A~~dz = T ' ,l=1 ' 2 '"" ~ 3 


Where, Au is the area of the up tube. 


6.2 The Reaction Kinetics 

The rate of reaction for the i* species depends upon the mass fraction of that species, 
mass fraction of oxygen and the Temperature field. The oxidation reaction of the i* 
species can be given as 
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A i + n0 2 — L - > pC0 2 + cjH 2 0 +.... 

Since this is an exothermic reaction, the heat of reaction would be negative. The continuity 
equation of the i* species is given by 


dco- m T dco- 

p — L +— i ; l - = r- 

dt A dz 1 


Where, A is area of up or down tube and m T is the local mass flow rate. The rate of 
reaction is given by 


N 


r i = ~ n i k o e RT n®f Pi 
1=1 


Therefore, the continuity equation becomes 


do): m T dco- 

» L j i L — 


N 


dt A dz 


~ ~ n i k o e RT Tlwf' 


or. 


dco- m T dco- 
p — l - +— i l - + 

dt A dz 


N 


n i k o e RT n what?' 


-l 


= 0 


or. 


dt A dz ' ' 


where. 


£ 

Q i =n i k 0 e RT fl wficofr 1 

M 


6.3 Variational Formulation 

The spatial descretization remains same as described in chapter 4. The Galerkin 
variational form can be used for the additional governing equations. The various 
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elemental matrices for a pseudo-element formed by one-dimensional real elements of up 
and down tube can be thus obtained as 
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The elemental force vector becomes. 
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f' 




/i+ b 2 ) f * 

Cl , Y2 c N 


+ - / 3 + t + r 


V3 2. 

1 c 
- + — 

3 2, 


f% + 


3 

A 

U 


/1 +|| + ^|/ 4 


The assembly procedure described in chapter 4 applies here too. 



Chapter 7 

Conclusion and Further Scope 
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7.1 Conclusion 

The Finite model uses Galerkin formulation for the heat conduction equation for 
earth and Petrov-Galerkin formulation for the convection equations of the reactor tubes. It 
uses the Predictor-Multicorrector algorithm for solving the time dependent ordinary 
differential equations for a more accurate solution. 

The model can be used to depict the time required for start up of the reactor and 
mass flow rate of the water required for a desired starting temperature. This model can 
also be used for comparing different geometries of reactors and efficiency of the reactor 
for different fluids. 

The results obtained here are in agreement with the results of the previous study 
made by Mittal S (1990), within the machine generated computational error, as different 
machines were used for simulation. 

7.2 Further Scope 

This study can be carried forward with making some generalization. The friction 
in the pipe has been neglected here. The friction can be included, and then the pressure no 
longer is hydrostatic. Also, the properties of fluid in the down and up tubes can be 
assumed to be varying with temperature and pressure. 
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